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TECHNICAL MEMORANDUM 


SPACE SHUTTLE MAIN ENGINE VIBRATION DATA BASE 

INTRODUCTION 


The Space Shuttle Main Engine (SSME) single engine test program has been 
underway since 1975. Throughout this time, engineers have been responsible for 
reporting on the dynamic characteristics of the engine and its components. The 
specific components of primary concern dynamically have been the two low pressure 
pumps and two high pressure pumps. Increased engine performance requirements, 
more aggressive testing schedules, and demand for higher resolution analysis have 
led to the development of the SSME Vibration Data Base. This data base has been 
operational since 1983 and has decreased data reduction turnaround time by a factor 
of three, while increasing significantly the diagnostic analysis capability. 

The primary objective of the SSME Vibration Data Base is to make available to 
the dynamicist, in an efficient, systematic and timely manner, data analysis techniques 
which can be used in the evaluation of the operational integrity of the SSME turbo- 
machinery. The vast majority of the evaluations can be performed with conventional 
power spectral density (PSD) and root-mean-square (RMS) time history analysis. For 
those cases where these types of analysis are not adequate or where failures or major 
incidents occur, the analog magnetic tapes are available for more detailed and elaborate 
analysis. Analog magnetic recordings of all dynamic measurements are obtained for 
each test as standard operating procedure. Further, during the history of testing of 
the Space Shuttle Main Engine (SSME), there has been a need to develop techniques 
to describe special (peculiar) characteristics of the SSME turbomachinery. As a 
result, special purpose diagnostic algorithms, which tend to be somewhat heuristic in 
nature, have been developed in order to define these characteristics and to provide 
insight into their fundamental mechanisms. 

Both the PSD /RMS analysis and the special purpose diagnostic algorithms are 
discussed in detail. These routines have evolved to the present state over a given 
period of time. Adjustments are continually being made in the algorithms in order to 
make them more efficient and accurate. As the need for more algorithms occur, they 
will be developed and incorporated into the data base. Collectively, these algorithms 
eliminate the routine drudgery associated with engineering data reduction analysis. 

This therefore allows the dynamicists more time to spend in assessing the effects of 
the dynamic characteristics of the components under test. 


DATA BASE COMPONENTS 


The SSME Vibration Data Base consists of three components: (1) Power Spectral 

Density (PSD) Isospectral Data Base, (2) Diagnostic Data Base, and (3) Anomalous Fre- 
quency Data Base. The Isopectral and Diagnostic Data Base will be described in 
detail in this memorandum. The Anomalous Frequency Data Base will be presented in 
a later document. Each data base component is a self-contained software system; 
however, the Isospectral Data Base has the capability of communicating data to the 
Diagnostic Data Base. 


Each data base contains information acquired from vibration accelerometers 
located on the High Pressure Fuel Turbopump (HPFTP) and the High Pressure Oxygen 
Turbopump (HPOTP). The Diagnostic Data Base contains data associated with sixteen 
measurements per single engine test or STS flight. The Isospectral Data Base can 
handle up to 75 measurements per test or flight. This extra set of data includes low 
pressure fuel and oxygen turbopump accelerometers, internal strain guages and 
accelerometers on the high pressure pumps, and other external measurements of 
interest. A typical list of measurements is shown in Figure 1. 

THFSE ARE the measurements AVAILABLE. 
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Figure 1. Typical set of measurements for a 
single engine test. 


PSD ISOSPECTRAL DATA BASE 


The Isospectral Data Base is housed on the UNIVAC 1100 series mainframe com- 
puter at Slidell Computer Complex, Slidell, Louisiana. This data base requires a 
large amount of mass storage, hence the need for a large computer system. The 
basic unit of this data base is the Power Spectral Density (PSD). Contiguous PSD's 
for the entire test duration are stored for each measurement. The PSD is a fre- 
quency domain representation of a time signal. The general frequency composition of 
a signal is described in terms of the spectral density of its mean squared value [ 1] . 
The PSD is defined as: 
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Gxx(f) = |X(f,T) | 2 


( 1 ) 


where X(f,T) is the discrete fourier transform of the time signal. There are several 
advantages of storing the data in the frequency domain PSD format rather than the 
time domain. The PSD format needs one-half the computer mass storage required by 
time domain. This is critical due to the high digitizing sample rate used. Also 
important is the fact that the root mean square (RMS) value of any frequency range 
can be easily computed by simply integrating the PSD over that range. The RMS 
value can be obtained by 


Grms = 



Gxx(f)df 


A significant disadvantage of PSD format is that all phase information is lost. How- 
ever, for routine single engine tests, all of the necessary analysis required for data 
review reporting can be derived from PSD formatted data. Since the data is in the 
frequency domain, any discrete frequency can be easily detected and classified in 
terms of its speed and amplitude. Specifically important is the classification of all 
discrete frequencies related to the pump rotor rotation speed. 

The Isospectral Data Base requires a considerable amount of mass storage per 
SSME test. The factors affecting this amount are the test duration and the number 
of measurements processed. For example, a 300 second test with 20 measurements 
processed would generate 16.4 million floating point numbers. Given a 500 sec test 
(i.e., typical STS flight duration) with 40 measurements processed, the amount of 
floating point numbers would increase to 53.2 million. Due to these mass storage 
demands, only three to five complete tests can be on-line simultaneously. More on- 
line tests can be made available by limiting the number of measurements on-line per 
test. An example of this would be to have on-line only those measurements associated 
with the HPOTP. This would reduce the number of measurements per test from 40 
down to 6 to 10. There is a tremendous amount of processing and display capability 
in the Isospectral Data Base. Each computational task is made more efficient due to 
the fact that the data is only digitized once. In the pre-data base processing system, 
the analog data would have to be converted to digital form everytime analysis was 
done, resulting in repetitious processing. Now, most processes or tasks can be 
executed in the amount of time required to simply digitize the data. This ability to 
store a completely digitized measurement is a significant factor in reducing the data 
reduction turnaround time. The following is a description of the various analysis and 
display capabilities of the data base software. 


PSD Ensemble Averaging 

An ensemble averaged PSD can be computed and displayed for any user defined 
time interval. The user has complete control over the plot bandwidth, the number of 
decades in the Y-axis, the start time of the averaging computation, and the number 
of averages to be used. Each PSD plot annotates the six highest peak values listing 
the frequency and amplitude. Also annotated is the composite RMS value (i.e., the 
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RMS value computed over the plot bandwidth) . The synchronous speed RMS com- 
ponent is noted when applicable. An example PSD plot is given in Figure 2. 


RMS Computation 

There are two methods used in the data base software to calculate RMS values. 
The first method simply integrates over a fixed frequency interval to compute what is 
termed a composite RMS. The second method computes RMS values associated with 
pump speed discrete frequencies. This method uses a sophisticated speed tracking 
algorithm to detect synchronous related harmonics and compute the associated RMS 
components. The speed tracking algorithm will be discussed in detail later. Both 
methods use the same display routine which allows the user both automatic scaling 
and manual scaling privileges. The user has the capability of plotting RMS values or 
pump speed frequency versus time (Figs. 3 and 4) and also RMS values versus pump 
speed (Fig. 5). 


Isospectral Display 

Probably the most dramatic plot available in the data base is the Isospectral Plot 
(ISOPLOT). This routine plots consecutive PSD's on top of each other in a "waterfall” 
form (Fig. 6). A hidden line algorithm is used to preserve the three-dimensional 
character of the plot. This particular display is useful in identifying discrete fre- 
quencies. In just a couple of pages, the signature of the data can be defined for the 
entire test. All pump speed related peaks as well as those not related to pump speed 
can be readily identified. The ISOPLOT display can also be skewed from left to right 
(Fig. 7) or from right to left (Fig. 8). This tilting of the display allows the user to 
identify relative changes in amplitude of a given discrete frequency. It also helps to 
identify closely spaced frequencies. 


Peak Identification 

This routine analyzes each PSD and documents the six highest peak frequencies 
and their associated amplitudes. The user has control over the start and stop times 
of printout, and the minimum and maximum frequency range for the peak search. 

This particular routine is useful in identifying closely spaced discrete frequency com- 
ponents, particularly those which tend to cross over one another. An example output 
is given in Figure 9. 


Anomalous Frequency Tracking 

As mentioned earlier, all pump speed related components can be identified. 
However, in some SSME tests, discrete frequencies are present which do not qualify 
as a pump speed integer harmonic. Those frequencies are termed anomalous fre- 
quencies. These anomalous frequencies can take the form of unknown "wandering 
frequencies" in which the frequency and its amplitude changes continuously during 
the test, or of the form of subharmonics related to subsynchronous whirl characteris- 
tics associated with instabilities of the HPOTP and HPFTP. Anomalous frequency 
tracking can be executed by simply rejecting the previously identified pump speed 
harmonics and retaining the next highest peak value. This peak value is searched 
for within a user defined frequency interval. When an anomalous frequency is 
present (i.e., its amplitude is above the background noise), the tracking routine 
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Figure 3. RMS time histo: 
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Figure 4. Component speed trace. 
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Figure 6. ISOPLOT (not skewed). 
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Figure 9. Peak identification. 



locks on to it and computes its frequency, RMS value, and the ratio of its frequency 
to the current pump frequency. The onset of the anomalous frequency can be easily 
seen by observing the frequency versus time plot. During the time in which the 
anomalous frequency is not present, or its amplitude is below the background noise, 
the plot has a great deal of variation, i.e., the tracker is simply picking up random 
noise peaks. At onset, however, the tracker locks on and the frequency versus 
time trace smooths out dramatically. An example anomalous frequency plot is given 
in Figure 10. 


Summarized Data Computation 

During each static firing or flight, it is necessary to define a single value 
(G RMS) to characterize each measurement response at each power level. These 
response values are used to define the operating environment for an individual pump 
during a particular test. These values, when taken together, will describe the over- 
all health of each pump subsequent to the following test. Previously, these values 
were determined manually by the dynamicist during each test evaluation exercise. 

This manual technique would be to define an average- maximum value sustained during 
each power level of the test duration from the RMS time history data. Considerable 
engineering judgement was required to insure that no spurious peaks or other short- 
lived responses were considered while the maximum sustained responses were identi- 
fied. These average- maximum values were then inputed into the Diagnostic Data Base 
by hand for use in statistical analyses or other analysis of the data. 

An algorithm was developed to simulate this engineering method by computing 
11-point floating averages over each power level duration. It has been observed that 
a sliding average of this order closely approximates the value that would normally be 
obtained by a dynamicist evaluating the test. The maximum average computed is pre- 
served for each measurement at all power levels and output in tabular form (Fig. 11). 
This provides a quick summary of the overall character of the test under evaluation. 
The exercise of ruling out spurious peaks or other short-lived responses is left to the 
engineer. A significant feature of this routine is the ability to transfer this sum- 
marized data to the MSFC located Diagnostic Data Base. This host-to-host automated 
transfer saves hours of manual data entry time and eliminates manual data entry error. 


Speed Histogram 

A recent development in the Isospectral Data Base is a routine which calculates 
and displays a histogram of the pump rotation speed. This plot shows the cumulative 
time in seconds that the pump speed was active. This display is useful in defining 
the amount of time spent at potentially detrimental speeds. The user has total con- 
trol over all plot scales and the frequency bin width. The user also has the option 
to include more than one test in the analysis. An example histogram is shown in 
Figure 12. 


DIAGNOSTIC DATA BASE 


The Diagnostic Data Base is located in-house at MSFC, ED24, on a Tektronix 
microcomputer. The data storage medium is a TransEra 40 Mbyte Winchester Disk. 
This data base contains the summarized RMS values computed in the above mentioned 
routine. There are currently 599 single engine tests, 4 flight readiness firings, and 
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Figure 11. Summarized data tabulation. 



Figure 12. Speed histogram 



25 STS flights available to the user. Sixteen measurements are maintained, eight for 
the HPFTP and eight for the HPOTP. Other test related parameters are housed, 
including engine, HFPTP and HPOTP serial numbers, power level profile, and pump 
speeds at each defined power level. The serial number information is also stored in 
a separate directory file, giving the user the capability to access test(s) by a par- 
ticular pump or engine identification. This option allows for direct comparison between 
different pumps or pump-build configurations. It also saves the user from the time 
consuming task of typing in all the desired test numbers. 

Several processing and display routines are available to the user. The more 
widely used functions are described below. It is important to note here that all pro- 
cessing can be done for the composite, synchronous, two * synchronous, three * syn- 
chronous or 4 * synchronous RMS amplitudes. 


Directory 

This option not only lists the tests available, but also lists the test date, the 
engine, HPFTP and HPOTP serial numbers, and the power levels for which data is 
available (Fig. 13). This on-line per test information printout is a detailed window as 
to what data are contained in the data base. 


Lifetime Power Level Activity List 

This function documents the complete lifetime hot firing activity for a user 
defined engine or pump. This lifetime is broken down and categorized by each 
individual test and further by each power level active during the test. Using this 
printout, the user can immediately identify the amount of time (in seconds) a par- 
ticular pump (including revisions, if desired) has been active for each power level 
and the total operational time. An example output is given in Figure 14. 


Summary Sheet Printout 

This option prints out the RMS values for a user defined set of tests in a 
block chart format. The chart includes the RMS values for all 16 measurements, 8 
per page (Figs. 15 and 16). This output provides the user with a tabular listing of 
the RMS amplitudes over any test range. This document is also used as the standard 
overhead viewchart for reporting pump vibration performance during engine data 
reviews and evaluations. 


Statistical Analysis 

Without a doubt the most powerful computation routine in the Diagnostic Data 
Base, the statistics option has been used extensively to develop statistical models 
characterizing normal and abnormal pump behavior and to establish and update engine 
cutoff redlines. The statistics software computes the mean, standard deviation, and 
the maximum RMS amplitude for all sixteen measurements. Also computed is the 
probability density /distribution function using a GRMS bin width of 0.5. The output 
of this package includes a tabular listing of mean, standard deviation, and maximum 
RMS amplitude broken down by test stand (Fig. 17), probability distribution (Fig. 18), 
and probability density plots (Fig. 19). For the distribution and density plots, the 
user can request an overlay of several classical statistical functions. The application 
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Figure 13. Directory output format. 
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Figure 14. Component lifetime analysis. 
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Figure 15. Summary sheet for HPOTP measurements. 
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Figure 16. Summary sheet for HPFTP measurements. 
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Figure 17. Statistical tabulation. 
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Figure 18. Probability distribution 
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Figure 19. Probability density 



of classical distributions is desirable to enhance data characterization since this per- 
mits continuous statistical definition and manipulation from discrete measurement 
observations. This routine has recently been enhanced to calculate the statistical 
parameters associated with the RMS level acquired by dividing synchronous by com- 
posite (Fig. 20). 


DATA ACQUISITION PROCESS 


The initial process in the SSME Vibration Data Base is the data reduction/ 
conversion from continuous analog form to digital PSD form. To accomplish this task, 
an automated data reduction system was designed utilizing the equipment and man- 
power at the Slidell Computer Complex (SCC) Data Reduction Laboratory. This site 
was chosen due to its proximity to the single engine test stands at National Space 
Technology Laboratory (NSTL). SCC has not only the necessary data reduction 
facility, but also has the mainframe computer facilities required to house the massive 
amount of data. 

The automated data reduction system has revolutionized single engine test and 
STS flight vibration data analysis. The previous data reduction system was located 
at MSFC and consisted of microprocessor based Fast Fourier Transform (FFT) analyzers 
and digital waveform analyzers. This equipment has no significant data storage cap- 
ability; therefore, analog data had to be redigitized every time any new analysis was 
done. Another hindrance to this system was the inability to transport the analog 
tapes from NSTL to MSFC immediately following a test. In fact, the tapes would not 
arrive at MSFC until approximately 9:00 a.m. the following day. The new system is 
designed so that data reduction and display takes place during the late night and 
early morning hours following a test. This means that by the time analog tapes would 
arrive at MSFC under the old system, the new system has provided the engineers 
responsible for data review a complete data package consisting of three to four times 
the normal output of the old system. 

The data processing flow for this system is defined as follows (Fig. 21): 

1. Data are recorded at the test stand on analog magnetic tape. Following the 
test, the analog tapes are transported by automobile to Slidell Computer Complex 
(SCC). 

2. The data are converted from analog to digital form and stored on disk in 
time domain. The sample rate is 5120 samples /second. It is important to note here 
that the data are contiguous, i.e., there are no gaps or missing samples within the 
data. 


3. The data are recalled from disk in blocks of 2048 points. This is equivalent 
to 0.4 sec of time data. 

4. Each block of 2048 points is then converted from time domain to frequency 
domain using the FFT and then converted to PSD form using equation (1). The 
factor of two is applied to account for negative fequencies. There is no window 
applied to the time data before conversion. From the above stated sample frequency, 
the resulting PSD bandwidth is 2.5 Hz. This yields a maximum frequency of 2560 Hz. 

5. Each PSD, 1024 points, is then stored to digital magnetic tape. Consecu- 
tive PSD's are written to tape for the entire test duration. 
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Figure 20. Sync RMS /composite RMS probability density. 




















(Steps 1 through 5 are performed in the data reduction room at SCC using a 
Hewlett Packard 5451C Fourier Analyzer. Six measurements can be processed during 
a single pass through these steps. For a 500- sec test duration, this process takes 
approximately 3 hr for every six measurement pass.) 

(The remaining steps take place on the UNIVAC 1100 computer.) 

6. The PSD digital magnetic tapes are rolled in to the UNIVAC and made ready 
for processing. 

7. Technicians at MSFC access the UNIVAC through satellite computer links 
and execute plot and analysis software from Tektronix graphic terminals. For a 
typical afternoon test firing, the first six measurements will be available for plotting 
at approximately 2:00 a.m. the following morning. 


SOFTWARE DESCRIPTION 


The following technical descriptions are provided for potential users of this 
software. All of the described code is written in ASCII FORTRAN on the UNIVAC 
1100 system at SCC. All of the algorithms described can operate on up to six mea- 
surements at a time. 


Synchronous Peak Detection 

Before any pump speed related parameters can be computed, the synchronous 
speed has to be defined. This is a critical step in any turbomachinery vibration 
analysis process, particularly for the SSME turbomachinery because of the wide 
operating speed range; i.e., 20,000 to 37,000 RPM. The periodic (discrete) com- 
ponents present in the vibration data are directly related to the turbopump's operating 
speed. Subroutine "INTPSD" accomplishes this task. The algorithm consists of the 
following steps: 

1. Search interval definition — MSFC engineers have developed equations using 
past history data that identify the approximate speed of a pump given the current 
power level (percent rated power level) . The search range is from the predicted 
speed minus 25 Hz to the predicted speed plus 25 Hz. This range is adjusted for 
HPOTP data to detect the fourth harmonic of the synchronous peak by simply multi- 
plying the limits by four (single engine test data only). Likewise, the low pressure 
oxidizer turbopump (LPOTP) is adjusted to detect the eighth harmonic. 

2. Spatial average computation — The PSD for each measurement is summed 
together to obtain a single, spatial averaged PSD. This averaging reduces the noise 
level and enhances the periodic components. 

3. Peak identification — Once the search interval is defined, the algorithm then 
locates the highest peak value within the search interval and saves its corresponding 
channel number. 
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Speed Estimation 


Due to a fixed delta frequency of 2.5 Hz, the calculation of pump speed using 
only the discrete frequency peak is accurate within 150 RPM. To obtain a higher 
degree of pump speed computation accuracy, a routine was developed which will pro- 
vide an estimate of the speed by using not only the discrete peak but its correspond- 
ing side bands as well. The side band information is necessary due to the leakage 
effect caused by the FFT. In other words, the pump speed will be exactly correct 
if and only if the discrete frequency is in the center of the analyzer spectral band- 
width . 


The interpolation equation used to provide a measure of leakage uses a coeffi- 
cient "R" defined as: 

4 


R = [(Ai - N)/(Aj - N)] , 


where : 


Ai = Maximum mean squared amplitude within speed search window 

Aj = Mean squared amplitude of left or right side band of Ai (whichever is 
greater) 

N = Background noise mean squared level. 


This coefficient is then substituted into an empirical equation which provides an 
estimate of the frequency of the discrete signal as a function of percent of analyzer 
bandwidth. This equation has been verified experimentally as well as theoretically 
and is used to compute the channel correction factor (CCF). This equation is defined 
by: 


CCF = (0. 1878855/R ** 3) - (0.2290493/SQRT(R ** 5)) + (0. 13210677/R **2) 

- (0.7820805/SQRT(R ** 3)) + (1.206314/R) - (0.01570345/SQRT(R)) . 


The CCF is added or subtracted, depending on which side band is used (upper or 
lower), to the peak channel number, then multiplied by the delta frequency to obtain 
a more accurate estimation of the true speed. 


Anomalous Frequency Detection 

Any dominant or lingering discrete frequency which is not a multiple of the 
pump speed is termed an anomalous frequency. As described, this routine tracks 
these frequencies by rejecting the synchronous speed harmonics of all engine compo- 
nents and analyzing the remainder. The routine was developed as a spinoff to the 
synchronous speed detection scheme. In fact, the speed detection routine is executed 
prior to the anomalous frequency search. Once the synchronous speed is determined, 
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the anomalous frequency routine searches for the maximum amplitude within the user 
defined frequency range. Any frequency point within two bandwidths either side of 
a synchronous harmonic is rejected. The recent development of the speed estimation 
routine has made it possible to include points closer to the synchronous harmonic due 
to the higher resolution speed computation. 

Once the maximum amplitude is found, its RMS value is computed using a user 
defined number of side bands as well as the peak. This will provide an accurate 
estimate of the RMS value and is required because of leakage effects. Its frequency 
is then determined, and the ratio of the synchronous frequency to the anomalous 
frequency is calculated. These three values are saved for later plotting. 


Standard Test Data Dump 

For each single engine test or STS flight, there is a standard data package 
which is output by the isospectral software. This output consists of: 

1. Selected PSD plots throughout test. 

2. Composite RMS time history. 

3. Synchronous RivIS time history. 

4. 2 * synchronous RMS time history. 

5. 3 * synchronous RMS time history. 

6. 4 * synchronous RMS time history. 

7. Synchronous speed time history. 

All of the above are executed for each measurement. Since these data are computed 
and displayed for every test, a routine was developed which calculates /displays all of 
the above automatically during one pass of the PSD data. Items 2 through 5 are cal- 
culated using two subroutines, "CSYNC" and "CRMS,” while the first item is deter- 
mined by interrogating the test power profile and automatically choosing desired plot 
slice times. At the end of the RMS plotting function, a summarized data computation 
is executed. Also executed for HPOTP data is a synchronous speed histogram plot 
routine . 


FUTURE DEVELOPMENT 


The main objective for single engine tests is to develop a system which will 
acquire all of the high frequency data real time. This system will be located at the 
test stand and will eliminate the time-consuming post-processing using six measure- 
ments per pass currently being done. Immediately after the test has been completed, 
the system will begin the task of recalling data and computing Fourier transforms. 
Once the entire set of data has been transformed, the Program Support Communica- 
tions Network (PSCN) will be utilized to transfer the data to MSFC and Canoga Park, 
CA. All of the data digitized during a test will be available for plotting /analysis 
approximately 3 to 4 hr after a test. Development of this system has begun at MSFC 
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using a MASSCOMP MC-5700 super microcomputer. This system utilizes separate 
processors for the CPU and data acquisition as well as an array processor. Further, 
the data will be retained in Fourier coefficient form (real /imaginary) rather than the 
current PSD format. This will greatly expand the diagnostic analysis techniques that 
the dynamicist can apply in the evaluation process. 


CONCLUSION 


The development of the SSME Vibration Data Base has greatly improved the 
quality (and quantity) of post-test vibration data reduction. The system has improved 
productivity immeasurably by digitizing and storing the entire test duration instead of 
very small time intervals. The development of sophisticated computation, tracking, 
and display software has provided engineers with deeper insight and more complete 
information as to the dynamic environment of engine components. 

The development of the Diagnostic Data Base has made test-to-test comparisons 
more efficient by providing efficient search mechanisms along with tabulated, plotted, 
and statistical outputs. The Isospectral Data Base has been a significant aid in the 
evaluation of time, amplitude, and frequency trends within a test. With the develop- 
ment of a real time data acquisition system both for single engine tests and STS 
flights, the SSME Vibration Data Base will continue to develop into an extremely 
efficient and powerful data processing, evaluation, and diagnostic tool. 
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